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Abstract 


Recently, gravity measurements were made on a tall, very stable televsio 
transmitting tower in order to detect a non-Newtonian gravitational force. 
This experiment required the upward continuation of gravity from the Earth s 
surface ^points as high as only 600 ra above ground. _ The upward 
continuation was based on a set of gravity anomalies in the vicinity of the 
tower whose data distribution exhibits essential circular symmetry a 
appropriate radial attenuation. Two methods were applied to perform the 

upward continuation - least-squares solution of a local harmonic expansion 


methods yield comparable results, and 
of 50 /tGal or better (1 /iGal = 10 8 
commensurate with the tower gravity 
accuracy of 20 /iGal), and enabled a 
gravity. As expected, such precise 


and least-squares collocation. Both 
have estimated accuracies on the order 
m/s 2 ). This order of accuracy is 
measurements (which have an estimated 

definitive detection of non-Newtonian 3 . ■ - . . 

upward continuations require very dense data near the tower. Less expected 
was the requirement of data (though sparse) up to 220 km away from the tower 
(in the case that only an ellipsoidal reference gravity is applied). 

1. INTRODUCTION 

The upward continuation may be summarized as follows. A set of surface 
gravity anomalies, circularly symmetric about the tower, serves as the set 
of boundary values. The upward conti nuation is based on Newtonian theory 
either through a local harmonic series expansion of the potential or through 
the use of least-squares collocation. At altitude, the (upward continued) 
GRS67 normal gravity and attraction of the atmospheric layer are added to 
the upward continued gravity anomaly. The result is the total gravity a 
would be observed in a strictly Newtonian world. A comparison with the 
gravity directly observed using a gravimeter offers a test of the underlying 

theory. 

2. THE MATHEMATICS OF UPWARD CONTINUATION 

If V denotes the earth's gravitational potential , then under Newtoni an 
potential theory, V satisfies Laplace's differential equation in free space. 

V 2 V = 0. ^ 

One solution to Laplace's equation (1) is a Fourier-Bessel series 
expansion in cylindrical coordinates, which is appropriate for a local 
representation of the potential. Thus since in the 

the gravity anomaly is also a harmonic function (it satisfies Laplace s 
equation in free space), it may be expressed as the following series (Morse 

and Feshbach, 1953, pp. 1259- 


- 1262 ) : 


Ag(/>,0,z) = XI 


C (k) J n (k/») 


e - kz - in 0 dk, 


( 2 ) 


where p is radial distance in the horizontal plane, 6 is azimuth, and z is 
altitude. J n is the Bessel function of the first kind and n th order, 
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the wavenumber in the radial direction, and the C n are coefficient functions 
to be determined from gravity anomaly data. Since the final evaluation of 

• 't are ^ riown « is along the vertical, p= 0, this problem is 
simplified by defining an azimuthal average: 


A g (/>.z) = ^ J Ag(^ 


r 00 

,6,z) de = j C 0 (k) 


J 0 (M) 



(3) 


This average coincides with the unaveraged Ag along the vertical: 


00 

Ag(0,z) = Ag(O,0, z) = f C n (k) e~ kz dk . 

•'o 


(4) 


The second equality in (4) follows by noting that at p = 0 all J n ’s in (2) are 
zero except Jq which is one. It remains, therefore, to determine only the 
function Cq from azimuthal ly averaged gravity anomaly data. One can 
determine at best a finite set of values of Cq from a discrete and finite 
data set. The integral (3), therefore, is truncated to some finite limit 
and discretized. The truncation is optimized if the discretization is in 
the form of a Fourier-Bessel expansion. After making the appropriate 
substitutions in replacing the integral (3) with the discrete summation (see 

Romaides et al_., 1988), the least-squares solution is obtained by the 
following: 


p=l 


M 

X 
1 


m 


A m J n 
m 0 


<Vn.> 


r-/C- 


exp 


( v 


>>) 


— > min., 


(5) 


where P is the number of points, M is the number of zeros, /c m are the zeros 

°| Jo, fp - Ag(£pR,0p,Zp) - Ag(R,z), = pj R, r is the maximum radius of 
surface data, z 0 , z p are the elevations of the tower base and surface 
points, and A m are the solution coefficients. 

The Fourier-Bessel upward continuation then proceeds in three steps. An 
initial expansion is done using 23 zeros and approximately 1800 anomalies 
out to 220 km from the tower. This expansion allows the resolution of 
half wavelengths down to about 5 km. A second expansion is then done on 
only the residuals inside of 5 km (from the tower) allowing resolution of 
half-wavelengths down to about 100 m. And finally a third expansion is done 
on the second set of residuals inside of 50 m from the tower resolvinq half- 
wavelengths down to about 20 m. Table 1 shows the results of the three 
steps, and Figure 1 is a contour map of the final set of residuals. 

The other upward continuation method used is least-squares collocation 
which is an optimal estimation method based on the validity of Laplace's 
equation. All important in LSC estimation is a good representation of the 
covariance function of Earth's anomalous gravity field. The covariance 
model employed for the tower experiment consists of the actual degree 
variances of the Rapp-1981 field up to degree 30 (1300 km wavelengths) plus 
a sum of 9 reciprocal distance models (see Jekeli, 1984) each covering a 
specific band of wavelengths shorter than 1300 km. These models are fitted 
to a sequence of periodograms based on the gravity anomaly data. Figure 2 
shows these periodograms along with the power spectral density of the final 
model. The model psd at wavelengths shorter than 300 m represents a rough 


135 


extrapolation. Figure 3 shows a plot of the collocation weights that are 
applied to ~270 data points; note the increase in the last two sets of 
weights which could be an artifact of oversampling. 

3. CONCLUSION 

The upward continuation of the surface gravity anomalies was done using 
two independent computation methods. The two methods are in excellent 
aareement (see Table 2) with both results clearly showing a departure from 
?he inverse-square law! The conclusion is that there is a dominantly 
attractive non-Newtonian component to gravity. Previous experiments had 
indicated a repulsive component to gravity. This inconsistency can e 
overcome by postulating the existence of two additional forces, one 
attractive and one repulsive. Our data do not contain adequate resolution 
to distinguish between the one force (scalar) and two force (seal 
models but are consistent with both. Figures 4 and 5 show plots of the two 
models along with our data. The error bars used are those of the first 

method. 
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Table 1. 


Table 2. 


FOURIER-BESSEL UPWARD CONTINUATION 


OBSERVED MINUS PREDICTED (/* Gal) 


ELEVATION STEP 1 


0.69 

.422 

7.58 

.256 

9.38 

.233 

23.07 

.094 

45.93 

-.054 

68.76 

-.159 

93.92 

-.187 

192.17 

-.341 

283.58 

-.455 

379.54 

-.534 

473.24 

-.567 

562.27 

-.566 


STEP 2 

STEP 

.266 

.009 

.116 

.002 

.098 

.008 

-.013 

-.013 

-.121 

-.098 

-.196 

-.176 

-.199 

-.184 

-.306 

-.301 

-.414 

-.412 

-.499 

-.498 

-.541 

-.540 

-.548 

-.547 


ELEVATION METH 1 ERROR METH 2 ERROR 


0.69 

9 

7.58 

2 

9.38 

8 

23.07 

-13 

45.93 

-98 

68.76 

-176 

93.92 

-184 

192.17 

-301 

283.58 

-412 

379.54 

-498 

473.24 

-540 

562.27 

-547 


59 

35 

9 

59 

-2 

13 

59 

0 

16 

58 

-27 

36 

57 

-100 

61 

56 

-171 

80 

54 

-179 

95 

48 

-304 

117 

44 

-413 

120 

38 

-493 

120 

37 

-528 

121 

36 

-526 

121 
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LOG* (GRAVITY ANOMALY PSD lMGAL/(CY/M|*) 
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